/ 


NASA TECHNICAL NOTE 


CO 

o 


lyO 


! 

; 

t 



ftf3- /S ^ xe> 

NASA TN D-7134 




ANALYSIS OF ELECTRICAL 
i CONTACT OCCURRENCES BETWEEN 
i ROLLING SURFACES WITH APPLICATION 
TO ELASTOHYDRODYNAMIC LUBRICATION 

• by Steven M. Sidik 

Lewis Research Center 
Cleveland, Ohio 44135 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION • WASHINGTON, 0. C. • JANUARY 1973 i 



2. Government Accession No. 


3. Recipient's Catalog No. 


1. Report No. 

NASA TN D-7134 


4. Title and Subtitle ANALYSIS OF ELECTRICAL CONTACT OCCUR- 
RENCES BETWEEN ROLLING SURFACES WITH APPLICATION 
TO ELASTOHYDRODYNAMIC LUBRICATION 

5. Report Date 

January 1973 

6. Performing Organization Code 

7. Author(s) 

8. Performing Organization Report No. 

Steven M. Sidik 

E-7166 


10. Work Unit No. 

9. Performing Organization Name and Address 

501-24 

Lewis Research Center 

11. Contract or Grant No. 

National Aeronautics and Space Administration 


Cleveland, Ohio 44135 

13. Type of Report and Period Covered 

12. Sponsoring Agency Name and Address 

Technical Note 

National Aeronautics and Space Administration 

14. Sponsoring Agency Code 

Washington, D.C. 20546 


15. Supplementary Notes 


16. Abstract 

A mathematical model enabling the relation of electrical contact occurrences between rolling 
surfaces to the mean thickness of an elastohydrodynamic lubricant film is discussed. A previous 
development by other authors is shown to contain an error invalidating their results. Some 
potentially useful approximations are developed for the case when the lubricant film is large 
with respect to the surface roughness. Complete solutions that are valid for all film thicknesses 
are not yet available. 


17. Key Words {Suggested by Author(s)) 


18. Distribution Statement 



Elastohydrodynamic lubrication 
Randon noise processes 
Rolling surfaces 

l 

Unclassified - i 

unlimited 


19. Security Classif. (of this report) 

20. Security Classif. (of this page) 

21. No. of Pages 

22. Price* 

Unclassified 

Unclassified 

25 

$3.00 


For sale by the National Technical Information Service, Springfield, Virginia 22151 
























ANALYSIS OF ELECTRICAL CONTACT OCCURRENCES BETWEN 
ROLLING SURFACES WITH APPLICATION TO 
ELASTOHYDRODYNAMIC LUBRICATION 
by Steven M. Sid Ik 
Lewis Research Center 

SUMMARY 

A mathematical model enabling the relation of electrical contact occurrences be- 
tween rolling surfaces to the mean thickness of an elastohydrodynamic lubricant film is 
discussed. A previous development by other authors is shown to contain an error in- 
validating their results. Some potentially useful approximations are developed for the 
case when the lubricant film is large with respect to the surface roughness. Complete 
solutions that are valid for all film thicknesses are not yet available. 


INTRODUCTION 

It has been shown in reference 1 that a mathematical model for the elastohydro- 
dynamic (EHD) lubrication of rolling surfaces can be facilitated by representing the sur- 
faces as two-dimensional, stationary and ergodic, Gaussian random processes. This 
model has been used to estimate the average lubricant film thickness from experimental 
measurements by electrical contact methods. This report points out a fundamental 
error in reference 1, derives the correct method of analysis, and briefly discusses the 
areas where further theoretical developments must be made to make the mathematical 
model more realistic. Some approximate methods are proposed for both one- and two- 
dimensional analyses when the film thickness is large compared to the bearing rough- 
ness. 

Besides the supporting data presented in reference 1, Williamson (ref. 2) has dis- 
cussed the applicability of Gaussian processes to the analysis of rough surfaces. 

A fundamental fault in several previous analyses of two-dimensional random sur- 
faces (refs. 1 and 3, for example) is the assumption that the properties of surfaces are 



directly described by one-dimensional profile traces. Nayak (ref. 4) has discussed the 
appropriate techniques for characterizing two-dimensional surfaces from one- 
dimensional profile traces under the assumption of statistically isotropic surfaces. He 
shows that a naive analysis assuming that profile statistics may be directly used is er- 
roneous. The bulk of Nayak' s results are applications or extensions of results obtained 
by Longuet-Higgins in two lengthy and excellent papers (refs. 5 and 6). 

In this report, we begin with the one- dimensional case to develop the correct 
method of analysis. Then a generalization of the electrical contact occurrence method 
is made to two dimensions. Although a complete solution is not developed herein, ap- 
proximate solutions are presented for both the one- and two-dimensional cases when the 
film thickness is large compared to the surface roughness. A brief discussion of the 
assumption of isotropic surfaces is also presented. 

It is appropriate to remark at this point that we consider only a somewhat simpli- 
fied version of the real problem. Some of the complications arising in real life are 
bearing spin, bearing sliding, viscosity effects, nonuniform lubricant thickness, 
temperature- related effects, and so forth. However, until an adequate solution to the 
simplified problem which ignores these effects is found, little progress will be made on 
the more complicated problems. 


ONE-DIMENSIONAL ANALYSIS 

Relation of Film Thickness and Surface Roughness 
to Electrical Contact Occurrences 

The basic problem under consideration here is that of describing the behavior of 
two rolling surfaces separated by a thin film of lubricant. The purpose of the lubricant 
is to separate the metallic surfaces in order to prevent bearing wear and thereby 
lengthen the effective life of the bearing. Under differing conditions of load, speed, 
surface roughness, and so forth, the effectiveness of the lubricant film in achieving this 
purpose differs. This report provides a means of experimentally determining the mean 
film thickness. A first step in determining the average film thickness from experiment 
data and analyzing its effectiveness was presented in reference 1. 

The situation is represented in figure 1. There are two balls, conveniently referred 
to as the upper ball and the lower ball. Under a condition of load, because of local 
elastic deformation of the ball surfaces, there is a region of contact where the surfaces 
can be considered as locally parallel. This region is called the region of Hertzian con- 
tact and is illustrated in further detail in figure 2(a). 
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The authors of reference 1 have shown that, within the Hertzian area, the two sur- 
faces can be closely approximated by two Gaussian random processes separated by an 
average distance h which represents the average lubricant film thickness. The upper 
surface is denoted y y (x). We use the notation y u (x) ~ N (p u , C u (r)) to denote that y u (x) 
is distributed as a Gaussian random process with mean p u and autocovariance function 
(acvf)C u (r). The lower surface is denoted by y^x) and we have y^(x) ~ N(p^, C^r)). 

We assume both y u and y^ to be stationary and ergodic. The distance between the 
surfaces at any point x is described by y u (x) - y^x). Thus, we consider the compo- 
site process defined by 

s(x) = y u (x) - y z (x) - h (1) 

This process is illustrated in figure 2(b). We have (ref. 5) 

s(x) ~N(0 ,C s (t)) 
o 

where C_(t) = C ,(r) + C 7 (t). Note that CJO) = o. Whenever s(x) > -h, the surfaces 
are not in the state of metallic contact. Whenever s(x) < -h, the surfaces are in the 
state of metallic contact. Obviously, the physical process cannot have s(x) < -h. We 
assume that s(x) < -h corresponds to a state of elastic deformation of the surface as- 
perities, and that the process of elastically deforming and reforming asperities has no 
appreciable effect on the characterization of the surfaces in contact by equation (1). 

In order to determine, experimentally, when the surfaces are in metallic contact, 
it is possible to employ the electrical contact method described in reference 1. We de- 
scribe this technique briefly as follows: 

With reference to figure 1, we might apply a voltage across the upper and lower 
balls. Thus, whenever the balls are in metallic contact, a current flow will occur 
through the circuit. A continuous recording of the voltage will then provide a record of 
whenever the balls are in metallic contact. It is important to note that no voltage will be 
recorded whenever there is metallic contact anywhere within the region of Hertzian con- 
tact. A voltage will be observed only whenever the lubricant film is unbroken throughout 
the region of Hertzian contact. 

We assume that the surfaces are rolling against each other at a constant velocity 
with no sliding or spinning. Then if a voltage- averaging device is used during a test, 
the time- averaged fraction of the applied voltage is a function of the distance- averaged 
fraction of rolled-over distance which corresponds to the state of metallic contact. We 
denote the distance- averaged fraction of electrical contact as T c< The no- contact 
fraction is thus given by 1 - T c . We wish to derive a relation between the average film 
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thickness h and the contact fraction T £ so that experimental measurements of average 
voltage may be used to estimate h. 


Derivation of Contact Fraction 


We begin by describing what happens during an electrical no- contact occurrence. 
Let the length of the Hertzian contact be denoted by d. For each excursion of s(x) 
above the level s(x) = -h, electrical contact is not broken until the left edge of the 
Hertzian contact includes the first point at which s(x) > -h. That is, electrical contact 
is not broken until the right edge of the Hertzian contact is d units of distance past the 
first point for which s(x) > -h. Electrical contact is reestablished as soon as the right 
edge of the Hertzian contact includes the first point for which s(x) = -h again. Thus, if 
the excursion above -h is of length L, there is a loss of electrical contact only for the 
length (L - d) + where (x) + = MAX {x, 0 } . 

Define the following random variables (rv) associated with the random process 


s(x): 


6 (h) 
a v ' 


N a (t,h) 


N a (t,h|6 

— d) 

N a (t,h|a 

>d) 

N(t, h) 



rv denoting lengths of excursions of s(x) above level -h 

rv denoting number of crossings of s(x) from below -h to above -h 
in interval (0, t) 

rv denoting number of crossings of s(x) from below -h to above -h 
in interval (0, t) where excursion is of length less than or equal to d 

rv denoting number of crossings from below -h to above -h where 
the excursion is of length greater than d 

rv denoting total number of crossings of level -h in either direction 


Let the random variables 6^(h), N^(t,h), and N^(t, h 1 5 < d) be defined similarly for ex- 
cursions and crossings from above -h to below -h. Cramer and Leadbetter (ref. 7, 
ch. 11) present rigorous proofs that these can, in fact, be treated as random variables. 
The following developments are presented in a heuristic manner but can be made rigor- 
ous through application of the results in reference 7: 

We denote the distribution functions of 6 a (h) and S^(h) by A b (6 
spectively. Although there has been much research directed toward identifying these 
distributions, they are not known except for special values of h and certain forms of 
the acvf C(t) (e.g. , see ref. 8). We assume the expectations of 6 a (h) and 6^(h) exist 
and denote them, respectively, as t and t^. 

Of primary interest are the two conditional expectations T^ and T^, defined by 


a ) and B h (6 b ), re- 
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T N =E[6 a (h)|6 a ( h )>d] 

V E [ 6 a< h) l 6 a (h)sd ] 

where d is the length of the Hertzian contact. Since 6 (h) is a positive rv, we have 

3 . 

*a ' Vr[6 a (h) > d ] + T c P r [ 5 a (h) S d ] 


For sufficiently large t, the following statements are true (« denotes approximately 
equal): 


N a (t,h)«N b (t,h) «lN(t,h) 

2 


N a (t,h|6 a W>d) 

N a (t,h) 



N a(t, h |6a {h)sd ) 

N a (t,h) 



Because of the stationarity and ergodicity assumptions, the approximate equalities be- 
come equality when t — Likewise, the following derivations involving approximate 
equalities become exact in the limit as t — °°. Thus, for sufficiently large t, the total 
no- contact distance is approximated by the number of excursions above the level -h 
which last longer than d, multiplied by the average amount by which such excursions 
exceed d. Thus, 


N 


1 - T. 






>d 


l)(T N - d) 


But we have 


t w N & (t, h)t & + N^tjh)^ 

“<‘a**b,^ 
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1 - T 


c 


N a (t,h|8 a (h)>d)( T N 


<‘a + *b> 


N(t,h) 

2 



(T N -d) N a (t,h|6 a (h)>d) 
l a + t b N a< 4 » h > 






Pr (6 a (h) > d} 


E { 6 a ( h)l 6 a (h) > d] - d 

E[6 a ( h )] + E[6 b (h)] 



( 2 ) 


In order to utilize equation (2), we must know the distribution of 6 (h) and the expecta- 

cL 

tion of S b (h). In general, these quantities are unknown. There are some asymptotic 
results for 5 b (h) as h - °°, however. Cramer and Leadbetter (ref. 7) also have a 
number of results concerning the expectation of the excursion length which are exact. It 
must be noted that these exact results are, unfortunately, inapplicable since the quantity 
of interest is E (L - d) + , not E(L). Before we proceed to consider the application of 
the asymptotic results, we present a comparison of the results obtained here and those 
given in reference 1. 


Discussion of Results of Reference 1 

Reference 1 and reference 9 were the first to discuss the relation of electrical con- 
tact measurements to the average film thickness. The mathematical and probabilistic 
arguments forming the basis for the discussions in reference 1 are given in reference 9 
and a series of subsequent progress reports (refs. 10 and 11). Their results for the no- 
contact fraction for the one-dimensional case are not in agreement with ours and some 
discussion of this disagreement is in order. 

In reference 9 there are actually four slightly different statements for the expres- 
sion for 1 - T . Each of these statements is given here and each is then discussed. 

The notation of reference 9 is changed to agree with the notation of this report. 

(1) Reference 9, page 13: '’Thus, the (time) average fraction of the no-contact time 
is equal to the probability of finding intervals longer than d where s(x) > -h. ” 
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(2) Reference 9, page 13: 1 - T c = Prjs(x) > -h|x Q < x < x Q + dj 

(3) Reference 9, pages 13 and 14: "The probability (given in item (2) above) is equal 
to the probability of finding intervals of length L > d within which the random process 
does not cross over the line s(x) = -h, whereas it crosses this line upward ahead of the 
interval and downward behind the interval. This probability can be computed if the dis- 
tribution of "spacings" between subsequent upward and downward crossings of the 
process over the line s(x) = -h is known. " 

(4) Reference 9, page 16: "From Pr jfi^h) > dj we get the probability of finding a 
spacing 6 (h) > d at an arbitrary point x by observing that this event requires 

(a) That s(x) > -h at. x 

(b) That, given (a), the x point be part of an interval of length L > d in which 
s(x) > -h at all points 

Thus, the average no- contact time fraction 1-T C is 1 - T c = |l - $(-h/o)Jpr|6 a (h) > dj. " 

Statement (1) can be easily proved false by a simple counter example. Consider a 
Gaussian process defined by 


s(x) = sin(— 1 + e(x) - h 
\3d/ 

where e(x) is a band-limited white noise of arbitrarily small variance. That is the 
s(x) process may be well represented by a sine wave of period 3d. This is illustrated 
in figure 3. Note that for each excursion above s(x) = -h, there is loss of electrical 
contact only for the distance (L - d) + = d/2. For each period of length 3d, therefore, 
there is a no- contact fraction 1/6, and hence 1 - T = 1/6. It is evident, however, that 
the probability that any x chosen at random lies within an excursion whose length ex- 
ceeds d is 1/2. But 1 - T does not equal 1/2, and hence the statement is false. 

Statement (2) appears to be correct if interpreted properly. If we let all x be 
equally likely (not clearly possible since it would involve an improper distribution), then 
consider the probability of choosing an Xq such that s(x) > -h for all xe(xQ, x^ + d). 
The requirement that s(x) > -h for all xc(xq, Xq + d) effectively subtracts the length d 
from the excursions above s(x) = -h and hence accomplishes the same thing as com- 
puting the conditional expectation E^6 a (h) ( 6 a (h) > dj. 

Statement (3) is almost identical to statement (2) and the same comments should ap- 
ply. 

Statement (4) is the one from which the actual evaluation of 1 - T £ is developed in 
reference 9. This statement suffers from the same inaccuracy as statement (1). 
Namely, there is no provision for subtracting d from the length of the excursion. This 
may be seen by comparing their result to ours. From arguments presented in refer- 
ence 9 it can be seen that 
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Computational Method as h Approaches °° 

One of the primary purposes behind studying the EHD process is to design bearing 
systems where there is little or no metallic contact between surfaces. This implies 
that primary concern be given to the processes where h is large compared to the sur- 
face roughness. That is, where h /a is large. 

s 

For this situation some results from Cramer and Leadbetter (ref. 7) and from Rice 
(ref. 12) and Kac and Slepian (ref. 13) are valuable. These authors consider (1) the 
distribution of the lengths between successive upcrossings of a high level and (2) the 
distribution of the lengths of excursions above a high level. Since the Gaussian process 
is symmetrical about the mean, it is evident that similar results apply to excursions 
and crossings below a very low level. 

With regard to item (1), let F £ (x) denote the distribution function of the interval be- 
tween successive downcrossings of the level s(x) = -h, and let 1/0 denote the mean of 

c 

this distribution. It is shown in reference 7 that 



= Pr(the interval between one downcrossing and the next of s(x)= -h is < x/0 c ) 
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Since 


e 


c 


c;’(o) 

2-tt 



we have 


6 


c 


0 


h — 00 


Hence, the expected length between downcrossings increases rapidly and the stream of 
downcrossings becomes Poisson in nature. 

With regard to item (2), let F e (x) denote the distribution function of the lengths of 
excursions below the level s(x) = -h and let d Q denote the mean of this distribution. 
Then it is shown in reference 7 that 


lim 
h - 


F e ( e e x ) = 1 " exp 



= Pr(excursion length is < 0 e x) 


Since for large h 


0 «- 



we have 


e 


e 



0 


Hence, the expected length of the excursion decreases rapidly. 

These two results may now be combined to give a limiting expression for T c . In 

particular, for sufficiently large h, we will have excursions below -h occurring quite 

infrequently; and each excursion causes an electrical contact whose duration will be, on 

the average, d + 0 . Thus, for sufficiently large h fa and t, 
c s 
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We note that Rice (ref. 14) has shown 


N b (t,h) 

~ 1 

’-c;’(o)‘ 

1/2 

exp 

1 

i 

s* 

NS 

1 

t 

1 *= 
i 

1 

U 

i 


_2Cs(°)J 


and hence 


1 

77 


■Cg(0) 


ni/2 


C s (0) 


exp 


-h‘ 


2C g (0) 




+ I 


277 


c;’(o) 


1/2 


o 

Since C_(0) = a , the variance of the Gaussian process, this equation may be rewritten 
s s 

and solved in terms of h /a . In order to simplify this, we assume 6 is small rela- 

S 6 

tive to d so that 



or 



(3) 


2 » 

In order to apply this equation, ct“ and C"(0) must be estimated. One method of 

O S 

estimating C''(0) is to note that this value is also the negative of the second spectral 
s 

moment. Thus, the spectrum could be estimated and its moments found. 

This completes the discussion of the one-dimensional analysis. The two-dimensional 
case is of much more benefit practically; so we now proceed to consider this, using the 
insights developed from this one-dimensional analysis. 



TWO-DIMENSIONAL ANALYSIS 


Derivation of Contact Fraction 

At this point we reiterate the observation that the Hertzian region is (for two balls) 
a circular area. The authors of reference 1 began their development with the assump- 
tion that the Hertzian region is an arbitrarily narrow and approximately rectangular 
area. They then developed what is essentially a one-dimensional theory. After con- 
sidering the one- dimensional theory, they then made a rather weak heuristic generaliza- 
tion to the two-dimensional theory. 

The results of Nayak (ref. 4) show that a naive analysis assuming that a profile 
trace in a single direction may be used directly is erroneous. He shows how the two- 
dimensional moments may be estimated from a one- dimensional trace under the assump- 
tion of statistical isotropy. Both Nayak (ref. 4) and Longuet-Higgins (refs. 5 and 6) 
consider such statistical descriptions as distribution of maxima and minima, distribu- 
tion of peak heights, mean summit curvature, and so forth. Longuet-Higgins also dis- 
cusses a method of determining a sequence of estimating functions which converge to 
the two-dimensional spectrum. Although these are important quantities in surface 
analyses, they do not directly contribute to the problem of electrical contact analyses. 

The statistical problem requiring solution before electrical contact analysis can be 
rigorously applied is now outlined. A complete solution to this problem is yet to be 
developed. 

The composite process of interest can be described similarly to equation (1) as 

s(x,z) = y u (x,z) - y ? (x,z) - h 


where 


y u (x,z) ~n(m u ,C u (t x ,t z )) 
y z (x,z) ~n(m z ,C(t x ,t z )) 

and hence 

s(x,z) ~ n(o,C s (t x , T z j) 
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where 


C s ( V T z> = C u< V V + C / T x> T z> 


Figure 4 presents a visualization of the metallic- contact areas in the plane parallel 
to the (x, z) plane at level -h. For simplicity we assume the x-axis corresponds to the 
direction of rolling. The area over which the Hertzian area passes as the surfaces roll 
over each other is bounded by the two horizontal dashed lines denoted z = Zj and 
z = Z 2 . The areas for which s(x, z) < -h are indicated by the enclosed dotted regions. 
For simplicity, we also assume that the process is statistically isotropic and that the 
Hertzian area is circular. These assumptions are discussed further in the section 
Isotropy Assumption. 

The left circle corresponds to the position of the Hertzian area when electrical con- 
tact is first broken. The right circle corresponds to the position of the Hertzian area 
when electrical contact is first made after passing over the "hill. " 

From figure 4 we see that the quantity whose distribution we need to know is X. 
Since the distribution of X appears too complicated to derive at the current state of 
knowledge concerning two-dimensional random processes, we will consider a potentially 

useful approximation for the case when h /a is large. 

s 


Computational Method as h Approaches °° 

As h/cr g gets large we might expect the same asymptotic Poisson character of the 
number of "valleys" to occur in the two-dimensional process as occurs in the one- 
dimensional process. In this case, the metallic- contact occurrences are more appro- 
priately illustrated by figure 5. This figure presents a visualization of the metallic- 
contact areas in the plane s(x, z) = -h. The areas for which s(x, z) < -h are indicated 
by the enclosed dotted regions. The area over which the Hertzian area passes as the 
surfaces roll over each other is bounded by the two horizontal dashed lines denoted 
z= Z l and z = Zg. The dashed circle on the left of the region denoted I corresponds 
to the position of the Hertzian area when electrical contact is first made. The point 
(Xm> z m ) denotes the coordinates of the point where contact is first made. The dashed 
circle to the right of the region I corresponds to the position of the Hertzian area when 
electrical contact is finally broken. The point (xg,z B ) denotes the coordinates of the 
point where contact is first broken. We let X denote the difference Xg - x M . If X is 
relatively small compared to the diameter of the Hertzian area, we may introduce the 
following approximation. It is reasonable to expect z^ and Zg to be quite highly and 
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positively correlated. Thus, assume z M and z B are perfectly correlated and that 
Zjyj is uniformly distributed over the vertical distance bounded by z = and z = Zg. 
In this case we have 


E(L) = 2E(L m ) + E(X) 
" 2E(L m ) 


If the shape of the Hertzian area and its dimensions are known, it is a simple matter to 
compute E(L m ). For example, if the Hertzian area is a circle of radius d/2, then 


E< l m> = 





where 


Z A = 


Z 1 + Z 2 


Some of the results of Longuet-Higgins (refs. 5 and 6) and Nayak (ref. 4) concerning 
the distribution of summit heights may now be employed if we assume the surface is 
isotropic. We first introduce some additional notation and results about two-dimensional 
random processes which will be needed. We again note that because the Gaussian 
process is symmetric, similar results apply to depths of valleys. For the sake of no- 
tational correspondence we will continue in terms of peaks and summits rather than 
nadirs. This is equivalent to discussing the process -s(x). 

We let ^p(k) denote the power spectral density function of a profile taken at an 
angle p along the (x,z) plane. Nayak shows that this power spectral density and its 
moments are all the information needed to describe the distribution of summit heights 
if the surface is isotropic. In particular, let m^ denote the i^ moment of g (k). 

That is, 

m i = j ^ g p (k)k 4 dk = c£ } (0) (yflf 
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where C^(0) denotes the i^ 1 


along p 
Let 


evaluated at zero. 

* = s/m J /2 = s/a £ 


s 


derivative of the autocovariance function for the profile 


a = 


m m. 
o 4 


m r 


Y = 


,.M-r 

[2(2a - 3)J 
[ 2(0 - l)(2o - 3)] 


and 


C 1 =— — 
1 20-3 


Let D gum be defined as the expected number of summits per unit area and 
p(s*) denote the probability density function of the summit heights in a normalized 
random process defined by s(x,z)/a . Then 

5 


sum 



and 


P(s*) 


= ^JL' c i (s * r 


3(2o - 3) 


|l/2 


t 3 tta/2 


,- 1 / 2 ( s *) 2 1 


2 ir 


a 


|i + *d3)J[(s*) 2 - ij 


a T /2 exp 

L3(a - 1)J 



*\ 

o(s*) 2 

[1 + $(y)] ► 
> 

v we 

2(a- 1) 
= z i " z : 


simply multiply D by z A to get the expected number of summits occurring within 
sixm n 
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the strip per unit distance along the x-axis, s^. We recall the summits are few and far 
between. Thus, the probability that two peaks occur within the Hertzian area is negli- 
gible. 

The average fraction of contact distance T £ is thus given by s h , the expected 
number of summits that exceed the level h(h/cr in the normalized process) per unit 
distance along x, multiplied by the average distance contact is maintained for each ex- 
cursion. That is, 


T c * s h 2E < L M> 


1 

3ttv^3 



p(s*)ds* 


This quantity may be numerically evaluated. 


Isotropy Assumption 

It is of interest to note that the assumption of isotropy may not be necessary. This 
would be of some importance in practice. Isotropy is quite reasonable when considering 
balls rolling on each other. However, balls rolling in a raceway will not, in general, 
exhibit isotropy because of the particular processes used to finish the surfaces. In par- 
ticular, raceway surfaces would typically be expected to have long grooves and hills 
running in the direction of grinding and much shorter grooves and hills in the normal 
direction. Assume that the grinding direction is parallel to the rolling direction of the 
ball in the raceway. If we assume the acvf is of the form 

C s (T x- T z )= g s exp ['“ T x + fjT z] ' 

then it is easily seen that the process defined by s(x,z/ /3) is isotropic. The only ef- 
fect such a transformation of variables has on the analysis is that the shape of the 
Hertzian area will be changed. In particular, a circle would be transformed to an 
ellipse. This produces no fundamental changes in the analysis of electrical contact oc- 
currences. 

An alternative method of considering nonisotropic surfaces which does not involve 
knowledge of /3 is to extend the results of Nayak to nonisotropic surfaces as outlined in 
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his paper (ref. 4). Such an extension does not involve any new principles. Some addi- 
tional work would need to be done regarding the expressions for D and p(s*). 

sum 


Prefiltering of Signals 

Whitehouse and Archard (ref. 3) have discussed briefly the presence of very low 
and very high frequency components of the spectrum and the resultant effect upon surface 
analyses. This deserves some further discussion in relation to the theory of electrical 
contact occurrences as developed herein. 

We first consider high frequencies. A typical profile that might result from a sur- 
face which contains primarily low-frequency components superimposed upon a signal of 
high frequency with small amplitude is illustrated in figure 6. It may easily be recog- 
nized that the high-frequency component contributes very heavily toward the number of 
"summits” contained in the signal. Not all these "summits" have physical significance 
with relation to electrical contacts, however. We may assume that since the asperities 
on the surface undergo elastic deformation, the small asperities will disappear when the 
larger asperities come into contact. Thus, the data recorded for description of the 
surface should be filtered so as to remove the high-frequency component. 

Figure 7 presents a typical profile that might result from a surface which contains 
primarily medium-frequency components superimposed on a signal with low frequency 
and low amplitude. In this case the low-frequency component will be reflected as very 
long and smooth "hills" and "valleys. " But if these are quite large with respect to the 
Hertzian area, it is more reasonable to assume that the bearings will roll up over the 
hills and down through the valleys rather than cause there to be very long contact and 
no- contact occurrences. Thus, the data recorded for description of the surface should 
also be filtered to remove such low-frequency components. 

Since there are little data available in sufficient detail, it is not clear that this 
problem need even arise in practice. The most serious problem is with the high- 
frequency signals, since they contribute the most to the number of summits. Whitehouse 
and Archard (ref. 3) present some evidence that at least some surfaces may be de- 
scribed as in figure 6. In addition, they also comment that in some instances the smaller 
peaks undergo plastic rather than elastic deformation and hence disappear rapidly with 
use. This lends support to the procedure of filtering out the high-frequency components 
and also presents a potential means for determining where the high-frequency cutoff 
point should be. 

In any event, what maybe considered "high" frequencies and "low" frequencies is 
open for research. Clearly, the cutoff frequencies for filtering must be related to such 
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things as the Hertzian area dimensions, the rolling speed of the surfaces, the load, and 
the elasticity of the materials. 


CONCLUDING REMARKS 

The purpose of this report is to show that reference 1 presents an erroneous de- 
velopment of the relation between average film thickness h and contact fraction and to 
present a correct development of the relation. This has been done. The general analy- 
sis that is valid for all values of h is seen to lead to an unsolved (so far) problem in 
random noise theory. The problem is that of describing the distribution of excursion 
lengths above a specified level. This has not yet been solved even in the one- 
dimensional case except for certain forms of the autocovariance function and certain 
levels. 

We have, however, presented approximate solutions for the case when the average 
film thickness is large. These approximations should prove of considerable value in 
practice since the area of primary interest in bearing analyses is the case when h/a g 
is large. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 20, 1972, 

501-24. 
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APPENDIX - SYMBOLS 


A h( 5 a< h >) 

B h( 6 b«) 

C(r x ), C(t x ,t z ) 

C'(r), C"(r),l 
C 1 (t) J 

D sum 

d 

E(A) 

E(A|B) 

F c’ F e 

h 

L 

m i 

N a (t,h), N b (t,h) 
N a (t,h|6 sd)\ 

N b (t,h|6 <d) J 

N a (t,h|6 >d),l 
N b (t,h|6 >d)J 

N(t,h) 

n(mj C(t)) 

Pr(A) 

p(s*) 

4 (k) 


distribution function of 5 a (h) 
distribution function of 5 b (h) 
auto covariance functions 
denote derivatives of C(t) 

average number of summits per unit area 
length of Hertzian area 
expectation of A 

conditional expectation of A given event B 

distribution functions of crossings and excursions, respectively 

average lubricant film thickness 

length of an excursion of random noise process above a given level 

xi_ 

i spectral moment 

random variable denoting number of crossings of s(x) from below 
(above) -h to above (below) -h in interval (0, t) 

random variable denoting number of crossings from below (above) 
-h to above (below) -h in interval (0,t) where excursion is of 
length less than or equal to d 

random variable denoting number of crossings from below (above) 
-h to above (below) -h in interval (0, t) where excursion is of 
length greater than d 

random variable denoting number of either-way crossings of level 
-h in interval (0,t) 

normal random noise process with mean level p and autocovari- 
ance function C(r) 

probability of event A 

probability density function of s* 

power spectral density of process corresponding to a profile taken 
at angle p 
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s(x), s(x,z) 
s* 

T N’ t c 
T c 

t a’ *b 
X 

y(x), y(x,z) 
6 a (h), 6 b (h) 
€(x) 

0 e 

M 

a 

$(x) 

Subscripts: 

a 

b 

B 

l 

M 

P 

s 

u 


random noise process in one and two dimensions, respectively 
normalized random noise process 

expected number of summits per unit distance in normalized random 
process whose heights exceed h 

conditional expectation of excursion lengths given 6 (h) as greater than 

a 

d and less than or equal to d, respectively 
contact fraction 

expected length of excursion above (below) level -h 

distance between center points of Hertzian area at make- contact and 
break- contact locations 

random noise processes 

random variable denoting length of excursion of s(x) above (below) -h 
white noise process 

mean of distribution of length between successive downcrossings of 
s(x) = -h 

mean of distribution of excursion lengths below s(x) = -h 

mean level of random noise process 

standard deviation of random process 

error function or cumulative normal distribution function 

crossings and excursions above -h 
crossings and excursions below -h 
break- contact occurrences 
lower ball 

make- contact occurrences 
profile statistics 

composite process s(x) or s(x,z) 
upper ball 
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positive part of x, that is, MAX (0,x) 
profile statistics in direction p 
"distributed as" 


"approximately equal" 
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Figure 1. - Idealization of lubricated balls in contact. 
(Scale of Hertzian contact greatly exaggerated. ) 


y u (x) - upper ball 

| Region of Hertzian contact d with surfaces locally parallel 





Rolling direction 


Figure 3. - Process six) - sin + e(x). 


□ No metallic contact 
E3 Metallic contact 
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